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Tj" I Effective Hamiltonian for the kicked double well system was derived using the Campbell- 

Baker-Hausdorff expansion formula. Asymmetric model for the kicked system was con- 
Q ■ structed. Analytical description of the quasienergy levels splittings for the low laying 

r) , doublets was given in the framework of the model. Numerical calculations confirm appli- 

cability of the proposed effective asymmetric approach for the double well system with 
the kick-type perturbation. 
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1 Introduction 

The connection between the semiclassical properties of chaotic systems and purely quantum 

^O ■ processes such as tunneling is a reach rapidly developing field of research nowadays. Our insight 

. ■ in some novel phenomena in this field was extended in the last decades. The most intriguing 

^ , among them are the chaos assisted tunneling (CAT) and the closely related coherent destruction 

of tunneling (CDT). 

The first one is an enhancement of tunneling in perturbed low-dimensional systems at 
relatively high external field strengths and high driving frequencies (in order the singlet- doublet 
^ ■ crossings to occur) [H El E]. This phenomenon takes place when levels of the regular doublet 
undergo an avoided crossing with the chaotic state |H[5]. At the semiclassical level of description 
one considers tunneling between KAM-tori embedded into the "chaotic sea". The region of 
chaotic motion affects tunneling rate because compared to direct tunneling between tori it is 
easier for the system to penetrate primarily into the chaotic region, to travel then along some 
classically allowed path and finally to tunnel onto another KAM-torus [HI E] . 

CDT phenomenon is a suppression of tunneling when values of amplitude and frequency of 
driving force belong to some one-dimensional manifold in the perturbation parameters' space [8]. 
This phenomenon occurs due to the exact crossing of two states with different symmetries from 
the tunneling doublet. In this parameter region tunneling time diverges which means the total 
localization of quantum state on the initial torus. 
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CAT phenomenon as well as CDT were experimentally observed in a number of real physical 
systems. The CAT observation between whispering gallery-type modes of microwave cavity 
having the form of the annular billiard was reported in the Ref. [H]. The same phenomenon for 
ultracold atoms was experimentally investigated in Refs. [IHl [11] . The study of the dielectric 
microcavities provided evidences for CAT as well [12]. Both phenomena were observed in two 
coupled optical waveguides [I3l [H]. Recently experimental evidence of coherent control of 
single particle tunneling in strongly driven double well potential was reported in Ref. [15j . 

The most common methods which are used to investigate the CAT are numerical methods 
based on the Floquet theory [6l UHl [17]. Among other approaches to CAT we would like to 
mention the scattering approach for billiard systems [181 HH] and quantum mechanical ampli- 
tudes in the complex configuration space [201 I2ll [22] • There is an analytical approach based 
on instanton technique, which was proposed in [221 IMl I2S] and independently used in PSJ- Al- 
ternative approach based on quantum instantons which are defined using an introduced notion 
of quantum action was suggested in [27] • 

We will investigate a quasienergy spectrum of paradigmatic model for different physical 
systems, namely double well potential. We consider this system with perturbation of the 
kicked type. One of the most attractive features of kicked systems is the well investigated 
simple quantum map which stroboscopically evolves the system from kick n to kick n + 1 and 
greatly facilitates theoretical analysis. 

The main idea of this paper is to study the possibility to construct an effective autonomous 
model for the non-autonomous perturbed system using Campbell-Baker-Hausdorff expansion 
formula and to test it in numerical calculations of the quasienergy spectrum. Both CAT and 
CDT are connected with the behavior of the quasienergy spectrum (avoided or exact crossing 
of the levels). Thus the development of the methods for calculation of this spectrum is impor- 
tant for extending one's knowledge in CAT and CDT. We regard the kick-type perturbation 
which is proportional to x. This perturbation, in contrast to perturbation proportional to x^, 
destroys the spacial symmetry in the system which is important for presence of the CAT phe- 
nomenon [2B]- The main role in quantum dynamics of our system is played by the classical 
asymmetry. There is no chaos induced processes in it but in our future work we will use this 
approach to system with CAT and CDT. 

In this paper we propose the effective model for the kicked double well system which gives 
a possibility to simplify the numerical calculations of the quasienergy spectrum and allows to 
determine both analytically and numerically the quasienergy splitting dependence on both the 
perturbation strength and frequency. 

2 Effective Hamiltonian for the kicked system 

Now lets construct the effective Hamiltonian for the double well system with the perturbation 
of the kick-type. Hamiltonian of the particle in the double-well potential can be written in the 
following form: 

Ho = - \- ao x^ - a2 x'^ , (1) 

Zm 

where m - mass of the particle, oq, 02 - parameters of the potential. 



We consider the perturbation of the kick-type which is proportional to x 
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Vper = ex y. S{t — nT), 



where e - perturbation strength, T - perturbation period, t - time. 
Full Hamiltonian of the system is the following: 

H = Hq + Vper- (3) 

Now we will construct an effective Hamiltonian for the system under investigation using the 
following definition: 

exp{—iHeffT) = exp{—iex)exp{—iHoT), (4) 

where RHS is a one-period evolution operator. We restrict our consideration by sufficiently 
small values of both the perturbation strength and period. Using the Campbell-Baker-Hausdorff 
expansion formula for the kicked dynamical systems we can rewrite the last expression (j4]) in 
the following way (2?j : 
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With the definition of g{z) formula (jS]) can be expanded in the following form: 
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where the expression under the integral is a map in the power k which for sufficiently small 
values of both the perturbation strength and period can be rewritten as follows: 
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Substituting (jTj) in expression ([6]) one obtains the following form of the effective Hamiltonian: 

Heff 



P , 4 2 . ^'^ . ^ I r^f 2 rri rr2\ 

Vaox - a2X + —x + —p + 0{e ,eT,T ). 

The fouth term in RHS of the last expression has the same order in perturbation parameters as 
first three main terms. The fifth term is proportional to small parameter, namely perturbation 
strength. We restrict our consideration by terms without small parameters and neglect all 
terms with order higher than zero in the perturbations parameters. As a result we have the 
following effective Hamiltonian for the kicked double well system: 
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Heff = -z \- a^x^ - a2X^ + —x. [p) 

This is the Hamiltonian for the asymmetric double well potential without perturbation. In 
contrast to the kicked system it is autonomous. In the next section we will consider the 
properties for this system and construct the effective asymmetric model for the kicked system. 
In section H] the correspondence between kicked and asymmetric effective double well systems 
will be tested numerically. 



3 Effective asymmetric model for the kicked system 

Hamiltonian of the classical particle in the asymmetric double-well potential is the following: 

^' ~ ^^ (9) 
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where oq and 02 are parameters of the potential, a - asymmetric parameter. The asymmetric 
potential with parameters m = 1, ao = 1/128, 02 = 1/4 and a = 0.15 is shown in the figure [1] 
(thick solid line). Eight lowest energy levels are shown in the same figure with thin solid lines. 
The form of the last term in the RHS of the Hamiltonian (Q is more handy due to special 
choice of the asymmetric parameter a. Parameter a for this form of the Hamiltonian is equal 
to shift between bottoms of the wells (dashed lines in the figured]). 
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Figure 1: Asymmetric double well potential (thick solid line) with eight lowest energy levels 
(thin solid lines). Minima shift {a = 0.15) is shown by dashed lines. The model parameters are 
m = 1, ao = 1/128, 02 = 1/4. 

In order to construct the effective model for the kicked system we should to introduce an 
effective value of the shift in the kicked system. Comparing obtained effective Hamiltonian (^ 



and Hamiltonian for the asymmetric double well system 
parameter in the following way: 
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Now we can give definition of the proposed effective model: asymmetric double well system 
with the effective parameter cXe// defining by expression flTUl) is an effective model for the kicked 
system with the perturbation parameters e and z/. The parameters oq, ^2 and m are the same 
for both systems. 

The definition of the effective asymmetric parameter (TTOI) shows that perturbation strength 
and frequency appears in it as a product eu. This is the first advantage of the effective approach. 
We have effectively only one asymmetric parameter (a) instead of two perturbation parameters 
(e and u). The second advantage is the more simple way of the numerical calculations which 
will be discussed in the next section. 



The third advantage of the proposed approach is that sphtting for the doublets laying 
below the barrier hump in the asymmetric double well can be described analytically. The 
asymmetric model can be considered as a pair of shifted harmonic oscillators. The shift is 
equal to asymmetric parameter a. It is obvious that in this constructed system non-degenerate 
energy doublets have splitting A = cr. In asymmetric system splitting between levels remains 
close to a for all doublets lying below the barrier top as for case shown in the figure [1] This 
correspondence can be used in order to give analytical description of the quasienergy spectrum 
in the kicked system. Using expression fIlOp we can write down the formula for low laying 
quasienergy doublets' splittings of the time- dependent system (jS]) 



A = ./f-. (11) 

It worth to mention the linear dependence of the levels splitting on both the perturbation 
strength and frequency. The applicability of this analytical description will be tested in the 
numerical calculations in the next section. 



4 Numerical calculations 

For the computational purposes it is convenient to choose the eigenvectors of harmonic oscil- 
lator as basis vectors. In this representation matrix elements of the Hamiltonian ([1]) and the 
perturbation (|2]) are real and symmetric. They have the following forms (n > m): 



tto — a 
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huj [n + -\ +- [-gaQ (2m^ + 2m + 1) - a'.^{2m + i; 



+ 5m+2 n - {g ao{2m + 3) - 4) A/(m + l)(m + 2) 

+ Sm+A n ^ V(m + l)(m + 2)(m + 3)(m + 4), 
7 



Xmn = Sm+ln \l ^ Vm+T, 

where g = h/mu and a'2 = a2 + mu'^/2, h is Planck constant which we put equal to 1, a; - 
frequency of harmonic oscillator which is arbitrary, and so may be adjusted to optimize the 
computation. We use the value u = 0.2 with parameters m = 1, ao = 1/128, 02 = 1/4. The 
matrix size is chosen to be equal to 200 x 200. Calculations with larger matrices give the same 
results. System of computer algebra Mathematica was used for numerical calculations. 

In order to obtain quasienergy levels {rjk) in the kicked double well system directly we cal- 
culate eigenvalues (A^) of the one-period evolution operator e~*^^e~'^ and express quasienergy 
levels through the definition rjk = i In Xk/T. Then we get ten levels with the lowest one-period 
average energy which is calculated using the formula {vi\ Hq + V/T \vi) {\vi) are the eigenvectors 
of the one-period evolution operator). The dependence of the quasienergies of these ten levels 
on the strength of the perturbation is shown in the figure [H Quasienergies of the two doublets 
with the minimal average energy (thick lines in the figure [2]) has a linear dependence on the 
strength of the perturbation in the considered parameter region. They are strongly influenced 
by the perturbation while some of the quasienergy states are not. 
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Figure 2: Quasienergy spectrum for the ten lowest average energy levels. All levels are numbered 
in order of the average energy values. Solid lines - quasienergy levels for the kicked system. 
Thick lines - two doublets with the minimal average energy. Empty squares - shifted energy 
levels of the asymmetric model. The model parameters are m = 1, qq = 1/128, 02 = 1/4 
and u = 0.5. 
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Figure 3: Quasienergy splitting as a function of the strength (a) and frequency (6) of the 
perturbation. Filled circles - results of the numerical calculations for the kicked system. Empty 
squares - shifted energies for the asymmetric model. The model parameters are m = 1, 
ao = 1/128, 02 = 1/4. 



To obtain the same levels in the framework of the effective model approach we have to 
calculate the eigenvalues of the Hamiltonian of the asymmetric double well potential (Q with 
the asymmetric parameter defined by the expression flTUl) . Then we shift obtained eigenvalues 
to zone (— f 5 1) in order to compare results with ones in the kicked system. The result of 
calculations on the base of this procedure is shown in the figure [2] by empty squares for two 
lowest doublets. Comparing results of direct and model calculations we make the conclusion 
that the levels of the doublets laying below the potential hump are correctly described by the 
effective model. 

Performed numerical calculations for the kicked and the effective system give the depen- 
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dence of the ground quasienergy splitting both on the strength (figJSt^a)) and the frequency 
(figJSl^b)) of the perturbation. Filled circles in figures correspond to kicked double well system, 
empty squares to numerical results obtained in the framework of the effective model. There is 
good agreement between splitting's dependencies for these two systems. They are linear as it 
was predicted by expression (ITT]) . It should be mentioned that all dependencies for the effective 
model was obtained from one series of the numerical calculations. We fix model parameters 
Co, 02, m and calculate numerically one set of numerical points for the dependence on the asym- 
metric parameter. In kicked system we should to perform one series of numerical calculations 
for every dependency. This is the first advantage which was discussed in the previous section. 

The second advantage which we should discuss after description of the used numerical 
methods is a more simple algorithm of the calculations. In the kicked system we should to 
calculate eigenvalues of the matrix exponents. This is more difficult task than in asymmetric 
model where we calculate the eigenvalues of the system Hamiltonian. 

Analytical result ( ITTj) which was put forward as third advantage of the method is plotted in 
the figures [3] (a) and |3] (b) by straight solid lines. Numerical points lie close to these lines. The 
agreement between numerical calculations and analytical expression (ITTj) is good (near 6%) in 
the parametric region considered. 

5 Conclusions 

Effective Hamiltonian for the kicked double well system was obtained using the Campbell- 
Baker-Hausdorff expansion formula. Effective autonomous asymmetric model for this system 
was constructed. This model is more convenient in numerical calculations than kicked one. 
Results of numerical calculations show that model correctly describes quasienergy spectrum of 
the kicked system for low laying levels. 

The analytical formula for the ground quasienergy splitting dependence on both the per- 
turbation strength and frequency was obtained in the framework of the effective asymmetric 
model. This formula predicts linear dependence of the ground quasienergy splitting on these 
parameters for the small perturbation strength and period values. Numerical results for the 
quasienergy splitting as a function of the perturbation frequency and strength demonstrate 
linear dependence as well. They are in a good agreement with the formula (TTTj) . Proposed 
approach will be used in future for investigation of the CAT and CDT phenomena. 
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